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Abstract 



Various deformations of the position-momentum algebras opera- 
pg \ tors have been proposed. Their implications for single systems like 

the hydrogen atom or the harmonic oscillator have been addressed. 
In this paper we investigate the consequences of some of these alge- 
bras for macroscopic systems. The key point of our analysis lies in 
0^ ■ the fact that the modification of the Heisenberg uncertainty relations 

present in these theories changes the volume of the elementary cell in 
^^ , the hamiltonian phase space and so the measure needed to compute 

partition functions. 

The thermodynamics of a non interacting gas are studied for two 
Qh. members of the Kempf-Mangano-Mann (K.M.M.) deformations. It is 

(— ! \ shown that the theory which exhibits a minimal uncertainty in length 

predicts a new behavior at high temperature while the one with a 
minimal uncertainty in momentum displays unusual features for huge 
volumes. In the second model negative pressures are obtained and 
mixing two different gases does not necessarily increase the entropy . 
This suggests a possible violation of the second law of thermodynam- 
ics. Potential consequences of these models in the evolution of the 
early universe are briefly discussed. 

Constructing the Einstein model of a solid for the q deformed 
oscillator, we find that the subset of eigenstates whose energies are 
bounded from above leads to a divergent partition function. 
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I Introduction 

Modern physics is an edifice in which every stone is tightly linked to the 
others. A slight modification in one area may produce important changes in 
different fields. Quantum mechanics starts with the commutation relations. 
Once they have been fixed, using a representation of the Heisenberg algebra 
we can in principle solve the Schrodinger equation whose solutions poten- 
tially contain all the physics of the ( non relativistic) system studied. The 
correspondence principle which states the link between the commutator of 
two variables of the phase space and their classical Poisson bracket is one 
of the basic axioms of quantum mechanics. It is not deduced from another 
assumption and cannot be judged in an isolated manner : only the whole 
theory can be confronted to experience through its predictions. 

Some authors have investigated the consequences of the alterations of 
these commutation relations on the observable of some physical systems |l| 



- [nj. In particular, some deformations of the canonical variable commu- 
tators studied by Kempf-Mangano-Mann(K.M.M) induce a minimal uncer- 
tainty in the position (or the momentum) in a very simple way, providing 
a toy model with manifest non locality. The implications of some of these 
quantum structures have been studied for the harmonic oscillator JT|, |3|] and 



the hydrogen atom fL2| . The transplanckian problem occurring in the usual 
description of the Hawking mechanism of black hole evaporation has also 
been addressed in this framework, for the Schwarzchild and the Banados- 
Teitelboim-Zanelli(B.T.Z) solutions || |(J. 

The purpose of this paper is the investigation of the modifications these 
models induce, not in the characteristics of a single particle, but in the behav- 
ior of a macroscopic system. The modern presentation of thermodynamics 
basically relies on statistical physics. According to the system under study 
(isolated, closed or open), one uses an ensemble (microcanonical , canonical 
or grand canonical) and the corresponding potential (entropy, free energy, 
grand potential) to derive thermodynamic quantities (pressure, specific heat, 
chemical potential, etc.). 

What has this to do with the commutation relations? The thermody- 
namic potentials are related to the derivatives of the partition function and 
the partition function itself is an average (of a quantity which depends on 
the ensemble used) on the phase space. To define a measure on the phase 
space, one needs to know the extension of the fundamentals cells. In the 
wrongly called "classical" statistical Physics, the Heisenberg uncertainty re- 



lations are used to show that this volume is the cube of the Planck constant 
and the indiscernability of particles justifies the Gibbs factor. The difference 
with quantum statistical physics lies in the fact that in the last case the sum 
giving the partition function has to be made on a discrete set of states. For 
example, to explain the specific heat of a solid, the Einstein model treats 
each atom as a quantized oscillator. 

If one modifies the commutators, one changes the energy spectrum and 
the Heisenberg uncertainty relations. The measure on the phase space is no 
more the same and this results in new partition functions and consequently 
different thermodynamic behaviors . 

The paper is organized as follows. In the second section we give a brief 
remind of the basics of statistical physics (to fix the notations) and their ap- 
plication in the derivation of the thermodynamic properties of some simple 
systems. The third section begins with a short overview of K.M.M. defor- 
mations of the momentum-position commutator. The non interacting gas 
is studied in two special cases: when a minimal uncertainty in position or 
in momentum is present , in the non relativistic and the ultra relativistic 
regimes. The Einstein model describing a crystal is also investigated. The 
last section introduces the harmonic oscillator of the q deformed algebra and 
uses it to work out the Einstein model of a solid in this particular theory. 

II Basics of statistical physics. 

Let us consider, in the usual theory, a system which is in contact with a 
large heat reservoir and doesn't exchange particles with the surroundings : 
it has to be studied in the canonical ensemble |L3], p~4]| . Its equilibrium state 
will be described by a fixed temperature and a fixed particle number while 
its energy will fluctuate around a mean value. Strictly speaking, for such a 
system , the particle number N is fixed once and for all. But, one knows 
that when phase transitions are not present, the descriptions given by the 
canonical and the grand canonical ensembles are very close. This will be used 
to compute the chemical potential in the canonical ensemble which is more 
tractable. Taking for example an assembly of non interacting non relativistic 
particles in a cubic box of length L, one solves the SchrSchro dinger equation 

ihd t ip = Hip (2.1) 



corresponding to the hamiltonian 

H=^- (2.2) 

2m 

with the boundary conditions ^(O) = if) {Li). The energy eigenfunctions are 
given by 

<fij;(r) = —j=exp(—iEt) exp(ik.r) (2.3) 

v V 



with 

2tt 
V = L 3 , fc = — (n s , n y , n z ) (2.4) 

the rij being integers. The dispersion relation reads 

The canonical partition function Z(T, V, N) is defined in term of the hamil- 
tonian operator H by 

Z(T,V,N)=Trexp(-^pj. (2.6) 

Using the complete set of states <j>%(f) displayed earlier, one finds 

z (r ^i) = X>p(-J^). (2.7) 

k 



The sum on k is performed on the discrete set given in Eq (|2.4j ). When the 
box is sufficiently extended and the temperature not too low, the preceding 
sum is very well approximated by an integral. Using Z = Z(T, V, N) = 
(1/N\)Z(T, V, 1) N one then obtains( with a trivial change of the variable of 
integration) 

Z{T,V,N) = J^J d 3N pd 3N qeM-E(q v , Pv )/kT) (2.8) 

with E(p,q) = p 2 /2m . The formula E(p v ,q v ) stands for the energy of the 
system when it is in a configuration in which the positions and the momenta 
of the particles are respectively p v ,q v . Finally 

V N 1 
Z= AM"A^ A=(ft 2 /27rmA;T)i (2.9) 



The Gibbs factor N\ is present only when the particles are undistinguishable. 
The free energy is related to the partition function by 

F(T, V, N) = -kT In Z. (2.10) 

In the variables T, V, N, the pressure P, the entropy S and the chemical 
potential fi are given by the relations 

dF dF dF 

while the internal energy is 

U = F + TS (2.12) 

and the constant volume specific heat reads 

Cv = %■ (2-13) 

For the case under studyone has 

, n . V [27rmkT\ 3/2 \ 



S = Nk 



5 , y f27imkT\ 3/2 ' 
2 +M \n{-J^) 



(2.15) 



C v = ^Nk (2.16) 

and the equation of state reads 

PV = NkT. (2.17) 

The ultrarelativistic gas obeys the dispersion relation E = c\p\ . The 
corresponding partition function 



87Tk 3 T , . v 



N\h 3N \ c 3 



V N . (2.18) 



leads to the equation of state given in Eq( j2.17| ) and the following expres- 
sions for the chemical potential, the entropy and the specific heat 



/' 



-fcTln 



S = Nki A + ln — 




C v = 3Nk. 



(2.19) 

(2.20) 
f2.21) 



The formula Eq( [2.8|) admits a semi classical interpretation. The system 
can classically be symbolized by a point evolving in the phase space. The 
probability for the system to be in a configuration in which the first particle 
is in the region g x ± Aqi,pi ± Api , the second particle in the region q 2 ± 
Aq 2 ,P2±Ap 2 , ■ • • is proportional to e( -E ®™ and proportional to the number 
of elementary cells contained in the volume of the aforementioned region. At 
the quantum level, the Heisenberg uncertainty relation ApiAqi > h/2 assigns 
to each elementary cell a volume h. The number of such cells contained in 
the region under consideration is 



n 



d 3 pjd 3 qi 

h 3 



(2.22) 



which justifies the h factor in Eq. Q2.8p . The Gibbs factor comes from the 
symmetrization( antisymmetrization) of wave functions necessary for the de- 
scription of a multiparticle bosonic(fermionic) state in quantum mechanics. 
At the semi classical level, it is introduced by hand to avoid an increase in 
entropy when mixing two formerly separated quantities of the same gas. We 
shall see in the next section that new commutators change the Heisenberg 
uncertainty relation and hence the measure needed to compute the partition 
function. 

We shall mostly interest ourselves to two specific cases in this work: the 
relativistic and the non relativistic non interacting particles. In these two 
cases,the hamiltonian H solely depends on the momenta. For all the systems 
we shall consider, we will be interested in the variations of the thermodynamic 
quantities induced by the uncertainty relations coming from the modified 
commutators. The new partition functions Z will be related to the usual 
ones Z* by a relation of the form 



Z^J 



N 



(2.23) 



J has to be computed for each modified commutators. The free energy then 
becomes 

F = F*-NkT\nJ. (2.24) 

The thermodynamic quantities are affected in the following way : 

1 (97 
P = P, + NkT 1 — (2.25) 

1 <9 7 

S = S* + NklnJ + NkT-— (2.26) 

J oT 

1 dJ 
fi = fi*-kT\nJ-NkT-j — . (2.27) 

.F*, -P*, • • • are the free energy, the pressure ... obtained from the usual theory, 
with the unmodified Heisenberg uncertainty. The new internal energy 

1 (97 

U = U, + NkT 2 1 — (2.28) 

leads to the following specific heat at constant volume 

IdJ 



C V = C* V + 2NkT-^ + NkT 2 



J_(dJ_Y I d 2 J 

~J2 [df + ldf 2 



(2.29) 



The last system we shall be interested in is a crystal. Two approaches are 
available : the Einstein and the Debye models. They give essentially the 
same prediction for the specific heat but the last one, which explains the 
properties of a solid by advocating an exchange of phonons between its sites, 
gives a better correlation with experiment near zero temperature. We shall 
nevertheless use the Einstein model because of its simplicity : each site is 
treated as a quantized harmonic oscillator vibrating at a frequency u. The 
partition function is computed exactly 

Z(T, V, N) = [2 smh(Uuj/2kT)}- N (2.30) 

and leads to the specific heat 

»-K9' k^S% (2 ' 31) 

which varies from zero to 3R for a mole as the temperature is increased. 
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Ill K.M.M. theory. 

The modified commutation relations proposed by Kempf, Mangano and Mann 
fTJ] are of the form 

[xi,pj] = iH6ij(l + ax 2 + ftp 2 ) (3.1) 

To be more specific, let us consider the case a = 0. The other commuta- 
tors read 

\Pi,Pj] = 

[xi,Xj] = 2ihf3(piXj - pjXi) (3.2) 

The last relation means that we have a "non commutative geometry" since 
translations along different directions do not commute anymore. The pres- 
ence of p 2 ,^ 2 in Eq. (|3.1| ) implies that the rotational symmetry, whose gen- 
erators are 

^ ij = 1 + pjp (*& ~ ^ ^^ 

is still preserved. 

One of the remarkable features of this theory is the existence of a minimal 
length uncertainty h\[P because of the new Heisenberg inequalities: 

A^A Pj > hSij f 1 + p JT((Ap k ) 2 + < p k > 2 ) J . (3.4) 

The momentum representation is found to be given by the operators 

Pi.ip{p) = Piip{p) (3.5) 

Xi .^(p) = ih(l + Pf)d p ^(p) (3.6) 

acting on a Hilbert space in which the scalar product is given by 



d 3 p 



< M > = / n , o^ FipWiP)- ( 3 - 7 ) 



The presence of a minimal length uncertainty implies that no position 
representation exists. The concept which proves to be the closest to it is the 
quasi-position representation in which the operators are non local : 

Xi = Zi+ y/pttw(-hy[pd u ) (3.8) 



y^ rT )tan(-z^ v /^%). (3.9) 



Pi -- 

When we refer to a system as a "relativistic" one , this will obviously 
be linked to its classical behavior. This precision is important because these 
deformations give rise to theories in which even if the the Poincarre group 
was a symmetry at the classical level, it is broken at the quantum level. 

The constants a,/3 which appear in Eq (|3.1|) are free parameters. What 
numerical values can they assume? Considering the case a 7^ it was sug- 



gested [15] that the minimal length uncertainty Ti^ffi should be of the order 
of the Planck scale. We will adopt a more liberal point of view here. The 
only constraint is that the deformed commutator must reproduce what is 
predicted by the usual theory and observed experimentally at our energy 
scale. Our attitude is inspired by recent works which have shown that a new 



physics may take place well before the Planck scale [16|]. As we shall see in 
the next sections, the behavior of a non interacting gas can be used to derive 
new, thermodynamic constraints on these parameters. Typically, the new 
theories will predict new equations of state or (and) new specific heats. The 
experimental values will be used to obtain bounds on the parameters of the 
model. 

When one introduces the parameters length and momentum scales a, (3 , 
it is straightforward that with the Boltzmann constant k , the light velocity 
c and the mass m of any particle, one can construct on purely dimensional 
grounds the characteristic temperatures l/(/3mk) and c/{kyf]3) ( the first 
is non relativistic, particle dependent while the second is relativistic and 
universal) and the critical volume a _3//2 . We are interested in what happens 
at these temperatures and volumes. 

III-A The case a = 

III-A.l The non interacting gas 

Let us study a non interacting non relativistic gas. Thinking in terms of 
quantum mechanics from the start, one has to solve the SchrSchrodinger 
equation for a particle in a box. This has to be done in the quasi position 
representation when (a = 0) because of the lack of a position representation 
. For simplicity, let us consider the one dimensional case i.e. Eq( |2.1| ): with 



p given in Eq. (|3.9|). One obtains the solution 



V>(t, -r e~ iEt exp I ±r^=5 arctan ^2m(5%E) j . (3.10) 

When the boundary condition ^>(£, £ = 0) = ip(t,£ — L) is imposed, one finds 
the dispersion relation 

En= 2^K tan {-^—) (3 ' n) 

( n being an integer) already obtained in [|1], |5|. This may lead to a cut off 
in order to avoid an infinite energy. However, this is not obligatory since 
the energy can become infinite only for very special values of the parameter 
/3.1f one does not have a cut-off in the momenta, the energy is no more an 
increasing function of the quantum number n. If the deformation parameter 
assumes a special value which forces one to consider such a cut-off, then once 
the number of particles N is fixed, the energy of the system is bounded. 

In the two cases, the integration on the momentum which is given below 
has an infinite upper bound. The one particle partition function 

CO 

Z(T, V, 1) = 5> _ * (3.12) 

n=0 

can be approximated by an integral if L is big enough. Introducing the 
integration variable p by p 2 = 2mE and replacing the length by an integral 
on position, we find 

Z(T,V,1) = ^Jdxdp Y ^- 2 e^ 2 l 2mkT . (3.13) 

One sees that in the new theory, one can keep the habitual dispersion re- 
lation and just modify the elementary cell volume. This could be anticipated 
with a semi classical reasoning. The Heisenberg uncertainty relation derived 
from the equation [x,p] = ih(l + (3p 2 ) reads 

AxAp> -(1 + Pp 2 ). (3.14) 

It assigns to the elementary cells of the phase space of the new theory a 
volume h(l+(3p 2 ) which replaces the usual value h. From this we find a simple 
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recipe when dealing with the semi classical approximation: it is obtained by 
keeping the classical dispersion relation but modifying the measure in a way 
consistent with the Heisenberg uncertainty relation. 

Let us now go back to the multi dimensional case. Due to the presence 
of three dimensional space dimensions, % is replaced by h 3 . The rotational 
symmetry results in the replacement of p by p 1 ( This comes from Eq. ( |2.1| ) 
)so that the new partition function reads 

III-A.2 The ultra relativistic non interacting gas. 

As it is shown in the appendix, the new partition function leads to a non rel- 
ativistic theory which differs very little from the undeformed case at habitual 
temperatures. This could be expected since the characteristic temperature 
T c is very high for a physically reasonable deformation parameter. The only 
sector where something new is likely to occur is at temperatures at least 
equal to T c . But, at those temperatures, particles are essentially relativistic. 
With the light velocity, the Boltzmann constant and the parameter /3, one 
can construct, on dimensional grounds, a temperature T cr = c/(ky/]3) . A 
new physics will be seen to emerge above this characteristic temperature. 
Working for simplicity in the ultrarelativistic limit, the expression of J is 
defined , similarly to the case treated in Appendix IV. A by the integral 

/•oo 

J= / Mx)dx (3.16) 

Jo 



with 

(j>{x) = 27T- 1 / W (l + P^x 2 ) (3.17) 



-3 



C 2 



with x = —p. 



c 1 

As can be seen from Eq.( |3.17| ), J has no volume dependence. This saves 



the relation pV = NkT. The presence of the temperature and the absence 
of the number of particles in the expression of J results in the fact that 
the entropy receives two contributing terms ( Eq. (|2.20|) ) while the chemical 
potential last contributing term vanish ( Eq.( [2.19| )). The internal energy 
remains unchanged is modified and by way of consequence the specific heat 
at constant volume . 
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The dominant contribution to the integral is located near the positive 
real x satisfying 



_ /3^-x 3 - 4/5^4-a: 2 - x + 2 = 0. 



(3.18) 



As explained in the appendix for the non relativistic case, a saddle point 
approximation is easily computed once when realizes that 



(f)"(x) = -2u~ 1/2 e-* [1 + 13 



2rn2 



k 2 T 



-x 



(3.19) 



The relevant zone for the saddle point approximation correspond to high 
temperatures. One then finds , when T » T cr 



(3.20) 



x = c/{U(3) 



Introducing the universal relativistic characteristic temperature T c = c/(kTy/]3), 
one finds 



J = e t, 
S - Nk I 1 + In 



1 



T 

T CT ) 2 9 /% 1 /4 

3 



/i 



-kT{\ + In 



N\ kT \T 



V f (87T) 1 / 3 hc\ fT ( 



N \ kT J V T 
and the equation of state takes the form 

It is easily found that an adiabatic process takes the form 

V = c te NT 6 
while the law of action of masses reads 



ya yb 

^■A^-B __ terp-6(a+b-c-d) 
vc vd ~ C 1 



(3.21) 
(3.22) 

(3.23) 



(3.24) 



(3.25) 



(3.26) 
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The dependence of the constants appearing in the last two equations on the 
characteristic temperature is obvious. 

The saddle point approximation used so far gives a picture which, even if 
it is qualitatively correct, is not the exact result. For example, when /3 — 0, 
Eq (|3.21| ) gives J = 0.85 while one knows the correct value to be J = 1. The 
exact result is computed numerically and plotted in the figures. From them 
one extracts the more accurate values of the <Tj and finds that o"o = 1 as it 
should. The presence of the V, N in the expression of S makes this theory safe 
concerning the mixing entropy and the the second law of thermodynamics, 
contrary to what happens in the case (5 = as we shall see in the next section 
. The equation of state takes the form 

p = 3p(l + g(T/T cr )) (3.27) 

The properties of the function g(x) are obtained numerically. In Fig.l is 
displayed the ratio of the new energy density and the habitual one, as a 
function of the variable T/T cr . 

Let us now give an estimate of the bound thermodynamics imposes on the 
deformation parameters. Consider the Helium whose specific heat at constant 
volume assumes the experimental values 12AJK~ 1 mole~ 1 (The undeformed 
theory assigns the value 12.47 to any non relativistic gas). With the help 
of the appendix, we can compute the specific heat in the new theory Cy = 
12.47(1 + a((3mkT)). The measured value tells us that 12.39 < C v < 12.59. 
At first order, cr(x) ~ x 2 . Assuming T = 300°K, this gives /3 > 10~ 45 . The 
unit in which beta is given is s 2 kg~ 2 m~ . One then finds a minimal length 
uncertainty 7 > 1CT 16 meters which is quite close to the bound derived 
from atomic physics considerations [|l^]. Working with another inert gas( 



the argon), we end up with the same result. This value of (3 can now be used 
to obtain a bound on the relativistic characteristic temperature, T cr > 10 9 i^. 
This is very far from the value T cr ~ 10 16 i^ which corresponds to a minimal 
uncertainty in position of the order of the Planck scale. Taking the modified 
relations used here as a "fundamental theory" i.e. not a phenomenological 
approximation( coming from string theory for example) , one sees that a new 
physics may appear well before the Planck scale. Its specific effects will be 
seen essentially in the early universe but astrophysics may also provide a field 
of application since the temperatures involved can be near the characteristic 
relativistic temperature whose lower bound is 10 9 K. The case /3 = has 
been seen to have a specific heat which is constant ( see -Eg.( |2.29j) ); the 
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argument we have used with the argon and the nitrogen to infer a bound on 
(3 in the case (3 /) , can not be repeated here to constrain a. One can use 
the departure of the equation Eq. ( |6.15| ) from the relation PV = NkT to find 
the allowed zone a < 10^ — 3) ( the unit is m~ 2 ) which is quite large . 

III-A.3 The Einstein model for a solid. 

When a = 0, the spectrum of a harmonic oscillator is found to display a 
quadratic term []J : 

E n = E Q + an + bn 2 



with 



B » = i fiw 47? + f + ii 



a = huj 



Ay/¥ 



16r 



(3.28) 

(3.29) 

(3.30) 
(3.31) 

' (2mM3)2 (3 ' 32) 

so that b vanishes when /3 — > 0. The free energy of a crystal with iV such 
oscillators at its sites is 



where 



Tiuj 
Ay/r 

1 



with 



F = N(E -kT\nK (T 7 (3)) 

an + bn 2 * 



n=0 



Its specific heat at constant volume reads 



C v = 3Nk 



IT 



K1JTJ3) T2 K 2 (T,{3) 
K {T,/3) K {T,P) 



kT 



-T 2 



(3.33) 
(3.34) 



' Ki(B,PY 



with 



K\ = -—— and K 2 



dT 



dT 



(3.35) 



(3.36) 
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A numerical treatment shows that the behavior of the specific heat in this 
model is very close to the one obtained in the undeformed case, for values of 
/3 in the range we found for the non interacting gas in the previous section. 
There appears to be a difference near the zero temperature but it is so tiny 
that it is experimentally undetectable . 

III-B The case (3 = 0. 

In the new context, the SchroSchrodinger equation in the position represen- 
tation leads to the energy spectrum 

En= i / m^ V 

2m yarctan ^aL J 

for non relativistic non interacting particles in a uni dimensional box of length 
L. This leads to a measure which could be inferred from the uncertainty 
relation 

AqAp> -{l + aq 2 ). (3.38) 

The partition function corresponding to the three dimensional case is 

1 f d 3N pd 3N q 

z= h^N\i (TT^F exp( : ~ E/ ' kT) - (3 ' 39) 

The measure of the phase space displayed in Eq. (|3.39|) lowers the influence 



of the large q i.e of the boundaries . 

III-B. 1 Non interacting non relativistic gas. 

In this case the quantity J defined in Eq. ( [2.23J ) is found to be the volume 
integral 

J= X -l %- . (3.40) 

VJ (l + aq 2 ) 3 K ' 



As can be seen from Eq.( 2.25 ), the fact that J has an explicit volume 
dependence breaks the relation pV = NkT. The last contributions to the 
entropy and the chemical potential vanish also because J (a, V) doesn't de- 
pend on the temperature and the particle number. The internal energy 
remains unchanged ( U* = U ) and by way of consequence the specific heat 
at constant volume is also unmodified. 
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Usually, the positions appear in the integrand of the partition function 
through interaction terms. The interaction between two particles depends on 
the vector separating them but this is not the case of Eq( |3.12| ). The situation 
may look similar to that of a perfect gas in a constant gravitational field since 
the particles do not interact among themselves but couple to an external 
field. In this case, the partition function is the integral of the quantity 
exp (— p 2 — mgz) /{kT). One sees that for the case , the factor (1 + aq 2 ) 
lacks the proper kT piece to allow such a parallelism. As a consequence, 
the thermodynamics generated by the K.M.M theory can not be interpreted 
as the usual one with a non trivial interaction. We will see that the same 
consideration applies to the q oscillator in the last section. 

The appearance of the norm of the position vector raises the question of 
the origin from which it is measured. This is related to the question of frame 
dependence which is not specific to the theories under study. In habitual 
statistical physics, the partition function involves the energy which is not 
the same for all observers, even galilean. The frame chosen is that of an 
observer at rest with respect to the container. In the case we are studying 
here, one has to specify a special point which will play the role of the origin. 
If for example we consider a spherical container, the center is a point with 
special status. 

The integral displayed in Eq( |3.4(J| ) is easily performed if we put the gas 



in a spherical container of radius R : 

? 1 

(3.41) 



Va 



R \ R 1 , , — . 

H -= arctan(vcm) 



1 + aR 2 ) 2 2 1 + aR 2 2^ 

The Eqs. , (|2.25|) , (|2.26|) , (|2. 2 71 ) show that the pressure, the entropy and the 
chemical potential are all affected by the new commutation relations. Let us, 
for simplicity, study two limiting cases : the small and the huge volumes. 

What happens when the volume can not be neglected in front of a~^ 3 ^ 2 ^? 
The expression of J becomes 

. 1 



J 



+ c 



3/2y a 5/2y5/3 



(3.42) 



o?l 2 V L a 

with 

a = 2.4674, b = 17.546, c = 78.6488. (3.43) 

Retaining the first two contributions, the pressure takes the form 
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This theory predicts an important change in the thermodynamic proper- 
ties for huge volumes; the correction to the pressure for small volumes given 
in Eq. ( |6.15| ) in the appendix is negligible. But the last expression shows 
that for very big containers, the pressure of a non relativistic gas may be- 
come negative. However, this does not contradict any experimental bound 
since for typical values of a like the ones of |T!J, [15| , one would need a 
recipient whose volume is greater than the earth . A negative pressure is a 
priori puzzling because it does not appear in everyday thermodynamics. But 
one knows it is an essential ingredient of the inflation scenario in cosmology 
fTT|| . We will come back to this point in the conclusions. 

Similarly, one finds for the entropy and the chemical potential in the case 
of important volumes : 

^=VM-,lM(-.)if^^V /2 ^ (3.1.) 




N \ h 2 ) 




» = -kT{ln{( — )±(^^J /2 }}. (3.46) 

The formula giving the entropy in big containers does not involve the 
volume itself and tells us that an adiabatic expansion is nothing else than 
a isotherm. The expression of the chemical potential shows that when the 
reaction 

aA + bB — > cC + dD (3.47) 

attains equilibrium, the following relation holds between the number of par- 
ticles Ni of the species i 

iV C iV £> _ T 3/2(c+d-a-b) ( o Ag) 

N%N b B K ■ > 

with t a quantity which does not depend on the temperature or the pres- 
sure. This is very different from the action of masses law in which only the 
concentrations are constrained by a similar relation. 

We now investigate entropy mixing. Consider two gases A, B contained 
in two compartments ( of respective volumes Va,Vb of a container separated 
by an impermeable wall). They are at the same temperature and pressure, 
their respective numbers of particles are Na, Ng .Now, let us remove the 
wall. In the non deformed theory (a = j3 = 0), the temperature and the 
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the pressure remain the same ( due to the conservation of energy and the 
equation of state). According to Eq.( |3.22|) , the increase of entropy for this 
process is non vanishing. 



In the theory we are studying, the equation of state specified by Eq.( |3.44| ) 



and the internal energy lead to the same pressure and temperature after 
mixing. However, the expression of the entropy is now given by Eq.( 3.45| ) 



For huge volumes, Eq. fl3.22|) gives that AS vanishes. This differs significantly 
from the non deformed theory. Let us consider two important but equal 
volumes Va = Vb = V and two equal number of particles Na = Nb = N of 
two different gases. The habitual theory states that the entropy mixing is 
2Nk\og2 while this K.M.M deformation predicts AS = 0. 

III-B.2 The ultra relativistic non interacting gas. 

J is found to be exactly the same as in the previous section so that compared 
to the undeformed case, the state equation, the entropy and the chemical 
potential change while the specific heat at constant volume remains the same. 



S = Nk { 4 + In 



or 



/ 3 hc 



^TXV^tlC 



3 / 2 ) N \ kT 



(3.49) 




,.= -*rin ^iR^ (3-50) 



The absolute value of the chemical potential is lowered compared to the 
habitual case. The formula Eq( |3.25[ ) is still valid, but with a different r. In 
summary, above the critical volume a~ 3//2 , the equation of state takes the 
form 

IY- = l + f(V/V cr ) p = 3p(l + f(V/V cr )) (3.51) 

for a relativistic gas; in the non relativistic case, the factor 3 in the last 
equation should be divided by 2. The function f(x) is plotted in Fig. 2 . It 
begins at zero for a vanishing x and quickly hugs to the horizontal asymptot 
y — — 1. In this model, the parameter a is fixed once and for all. If one 
considers an increasing volume for a fixed quantity of gas( fixed N) and a 
fixed temperature( fixed T), the last two equations show that p/p and p/T 
go to zero. 

An important issue for gases presenting a non trivial equation of state 
is the possibility of phase transitions. For example, a Van der Waals gas 
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exhibits a condensation plateau when it is compressed at (fixed temperature 

)under a critical volume ( obtained by Maxwell construction). In the phase 

diagram , this is linked to the existence of a critical point where the two first 

derivatives of the pressure with respect to the volume ( on the appropriate 

isothermal curve) are vanishing. From Eq. (|3.30|) one can isolate P and see 

that 

dP NkT .„,„,, d 2 P NkT 



m \ * fl{V/Vcr) W> = -v^ h{yjVer) (3 - 52) 

h(x) = -l-f(x)+xf(x) f 2 {x) = 2 + 2f(x)-2xf\x)+x 2 f"(x) (3.53) 

The functions fi(x) and f2(x) are plotted in Fig.3 an Fig.4 respectively. They 
show that there is no such a phase transition here. 

What we have done for gases is similar to Maxwell Boltzmann statis- 
tics because we did not symmetrize the wave functions for the bosons nor 
antisymmetrize for the fermions. At the fundamental level, one can ask if 
the concept of spin is relevant in these theories and, in the case the answer 
is positive, one may wonder about the spin-statistics theorem in the new 
context. 

The commutation relations given in Eq. ( |3.1|) and Eq. (|3.2| ) are not in- 
variant under a Lorentz transformation. As the spin of a particle is defined 
through the behavior of its wave function under this group, one apparently 
looses the notion of spin here. However, this is not obligatory true. For 
example, taking the ultimate structure of space-time to be a particular com- 
mutative geometry, the relevant group of isometries is not the Poincarre one 
but a q deformation of it. There is a new algebra but a spin has been defined 
in these theories and the wave functions for particles of spin 1/2, 1, etc... have 
been found [pl|. Although the question has not been addressed in K.M.M. 
theory, we expect a similar situation to occur. Considering as an illustration 
the replacement of the Dirac equation by 

(7% + m)V> = (3.54) 

with p^ given by the generalization in three dimensions of Eq. (|3.9|), the 



relevant symmetries of this equation should be linked to an inhomogeneous 
algebra since the operator is non-linear. We expect a generalization of the 
notion of spin which may conserve the spin-statistic theorem. This is not 
really crucial here since we are performing a Boltzmann like construction 
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which is known to be valid at high temperature and does not discriminate 
between bosons and fermions. What we called " relativistic systems" in 
the deformed theories are nothing than systems described by the habitual 
relativistic equations, but with the momenta replaced by their quasi position 
expression in the model where a vanishes, for example. This can in principle 
be made rigorous using the functional integral introduced in [[J since in this 
formalism one does not need to introduce neither the commutation relation 
nor the scalar product from the start. 

The theories studied here are unusual in many respects, we have seen 
that a minimal length uncertainty can exist for specific parameters. The fact 
that the commutator of two quantum operators is not simply proportional 
to the Poisson bracket of their classical counterparts creates a situation in 
which a quantity conserved at the classical level may not be conserved at 
the quantum level. This implies that one has to be careful when studying 
ensembles more complicated than the canonical one. 

We studied a system which exchanges only energy with the surroundings. 
When the number of particles is not fixed, one has to work in the grand 
canonical ensemble. However, it is known in usual thermodynamics that 
the two descriptions yield essentially the same predictions( equation of state, 
specific heat, • • •). A real difference appears essentially when describing phase 
transitions. 

Our treatment concerned only non interacting gases but we know they 
only approximate the real ones when the density is small. We think it is 
physically irrelevant to include interactions for habitual gases like hydrogen, 
ammonia, nitrogen, etc ■ ■ • . Taking for example the theory where only /3 ^ 0, 
we saw that a new behavior sets in at very high temperature. But at this 
energy scale, nitrogen does not exist anymore and the use of an interaction of 
the Lennar- Jones type is not justified; at very high temperatures, the existing 
species are elementary particles like photons, neutrino, leptons, quarks • • •. 
One may expect to have a different situation for the model in which only a 
is non vanishing ; one may conceive that the expansion of the universe may 
provide the necessary huge volumes. However, as the universe gets bigger its 
density decreases and the mean distance between the particles becomes too 
important ;they nearly don't feel the interaction. 
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IV The q deformed oscillator. 

A space with non commutative coordinates is called a quantum plane; its co- 
variant differential calculus has been investigated [jl|. A quantum group can 
be understood as a set of linear transformations acting on it. The commu- 
tation relations compatible with this structure |20| are quickly summarized 
below. 

The momentum and position operators P, X satisfy the commutation 
relation 

q 1/2 XP - q- 1/2 PX = iU. (4.1) 

The deformation parameter q is strictly greater than one; the unitary element 
U has been introduced to allow hermitian momentum and position : 



p+ = p x + = x u + = u-\ 

It is also subject to the conditions 

UP = qPU , UX = q- x XU. 

Introducing the extra parameter M , one defines the operators 

a+ = aU 2M + PPU M 

a = aU 2M + [3U- M P 
and finds that they obey a modified Heisenberg algebra 

aa + — q~ 2RI a + a = 1. 



(4.2) 

(4.3) 

(4.4) 
(4.5) 

(4.6) 



We first investigate what changes in the Einstein model of a solid . In this 



theory, the harmonic oscillator has been defined [0 as having the hamilto- 
nian 

H = cua + a (4.7) 

with a + , a the operators defined in Eqs.( |4.4|) , (|4.5|) and whose algebra is given 
by Eq. f3). 

A special role is played by the eigenvalue 



E n 



to 



-2M 



(4.8) 



of the hamiltonian : 
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a) There is an infinite subset of the spectrum which is bounded from above 

by E [oo] : 



1 _ a -2Mm 



K = u ~ „_ 2M ^ E M rn = 0,l,2,---. (4.9) 



For these states, a raises the energy while a + lowers it. 

b) The remaining part of the spectrum 

E n = u( \ M +(3h 2n \ n = -oo,...,+oo (4.10) 

(with (3 depending on the mass m and the frequency to) faces the re- 
verse, common situation with a + playing the role of a creation operator. 

c) For E = E[oo], a and a + do not change the eigenvalue. 
To give an example, the operators 

P = -i^-+ i l 7 , X = x (4.11) 

dx ^l _ q -2M 

U = q -?t* p + p *) (4.12) 

are found to obey the relations given in Eqs. ( |4.1| -4.3) which define the q 
deformed phase space. Introducing the operators a, a + as specified by Eq. 
( |4.4| - |4.5| ) and expanding the hamiltonian ua + a, one finds it consists of two 
terms. The first part is the hamiltonian of a free oscillator and the remaining 
part can be seen as an interacting term as a series in yh with h = e q — 1. 
The link between the parameters a, (3, 7 on one hand and the mass, the 
frequency of the unperturbed oscillator on the other hand are 



a = ±^^= 7 = ±V2muj 8= . (4.13) 

V2Mh V2muj ' 

It should be emphasized that the K.M.M. oscillator presents a spectrum 
which reduces to the habitual one as the deformation parameter is sent to 
zero. This is evident from Eqs. ( |3.29| ) and ( ^.30| ). The situation of the 



q algebra is different and so one expects it should lead to an entirely new 
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thermodynamics, even for a small value of the deformation parameter. Let 
us study a crystal whose sites are occupied by q oscillators. 
Constructing the partition function 






z(t, v,i)=J2 e ~^ + E e_lf ( 4 - 14 ) 



m=0 



one finds that the first sum diverges since E^ < E^]- How can we handle 
this? The extreme reaction would be to disregard the all theory as unphys- 
ical. Another possibility may be a renormalization of the partition function 
but this is questionable since we are not doing quantum field theory. An- 
other point of view is that these states may be interpretated as spurious and 
eliminated. We do not have a clear option between these possibilities for the 



moment. Introducing a = \lu(3(3q one rewrites the remaining part as 

Z(T,V,N) = e TT^W J^ e~^. (4.15) 

If a < 1, one has 

lim exp(-cr n ) = 1 (4.16) 

so that Z(T, V, 1) still diverges. The last option would lead us to drop also 
the states corresponding to high positive values of n. The question would 
then be to fix the index where the truncature should take place. 

If we discard the E^, we were left with E^] as the lowest state. As 
it corresponds to E n for n = 0, we think it should be kept in the physical 
spectrum. So, our truncature begins at n — 1 in this case : 

Z(T,V,l) = e TT^FW £ . e ~w (4.17) 

n=— oo 

When a > 1, the situation is reversed : 

Z(T, V, 1) = e v-rtxT Y, e ~ avF - ( 4 - 18 ) 

?1=0 

Writing q — 1 + h and retaining the first two terms of the binomial formula 
to approximate q n , one deals with a partition function similar to the one 
used in the K.M.M model. The same results apply, provided one chooses 
a value of h in the range which gives to f3 the values compatible with the 
bound given in section 3. Discarding the extra states has one nice feature: 
the thermodynamics is smooth in the deformation parameter. 
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V Conclusions 

Our investigation of the thermodynamics implied by two special K.M.M. 
deformations of the commutation relations have shown that a new physics 
appears at high temperatures in one case and for huge volumes in the other 
one. In the two cases, the universe should provide the best laboratory since 
the modified commutators can produce specific, visible signatures . In fact, 
the evolution of the universe will be affected in more than one aspect. Of 
course, the new equation of state of the radiation gas which fills the early 
universe influences the evolution of the scale factor through Einstein equa- 
tions. In addition, the new expressions of the entropy change the form of the 
adiabatic processes. As in standard cosmology the expansion of the universe 
is thought to be isentropic, this will enter into play. The fact that chemical 
potentials receive contribution corrections will affect the densities at equi- 
librium . Finally, the theory exhibiting a minimal length uncertainty may 
forbid the singularity present in the standard big bang scenario. A similar 
reasoning was advocated to argue that the Hawking evaporation of a black 
hole may halt without using the complementarity hypothesis ||. 



Let us say a few words about the negative pressures found in section|375 . 
In the inflation scenario, one uses a scalar field whose classical time evolution 
generates such a non positive pressure for a given period, because of a well 



chosen potential flTj • Eq (|3.21|) suggests that one may achieve such pressures 



with ordinary matter, provided the commutation relations are modified. But 
to achieve inflation, one needs the relation p < — p/3 and this not obvious 
from Eq.( |3.21| ). One also needs a mechanism to return to normal matter i.e 
positive pressures. Certain mechanisms of inflation are inefficient because 
their period of exponential growth never ends. This is the case in topological 
inflation where the core of a monopole, a cosmic string or a domain wall 
expands endlessly. It should also be stressed that here the negative pres- 
sures appear only when the volume gets big enough, contrary to standard 
inflation where it happens very early in the course of the universe. But the 
true challenge here, which may rule out all the model corresponding to the 
case j3 = 0, a 7^ is that there is no evidence for a preferred origin. A 
straightforward possibility is that this theory may be interpreted as impos- 
ing a maximal physical volume since going above it one gets into trouble. 
One of the most important novelties is the behavior of the entropy mixing in 
the a = case which seems to spoil the interpretation of the entropy derived 
from Eq. ( |2.11| ) as a variable counting the number of states. The breaking 
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of homogeneity makes it likely that one should consider Tolman (position 
dependent) temperatures like in curved space times. But in such spaces the 
fact that the quantity Ty/—g 0Q is constant is a vital information which seems 
to have no counterpart in the theory a^O. This model is also useless in 
cosmology since the existence of an origin is in contradiction with the homo- 
geneity hypothesis. One may use all the undesirable consequences of this the 
model to argue it is unphysical, like the quantum theory on a Moyal plane 
modifying the time-position commutator [x , x»] has been discarded because 
it leads to a non unitary theory. 

We plan to analyze the cosmological incidences of the model a = in 
a near future. The influence of the new physics in the early universe is at 
the center of a recent work using Unruh and Jacobson modified dispersion 



relations pl| . Our approach will differ by the fact the dispersion relation will 
come from an assumption on the quantum space time . The commutation 
relations studied here can be interpreted as phenomenological consequences 



of string or M theory [22], |23|. This paper suggests that string cosmology is 
not uniquely characterized by the evolution of the fields which appear at the 
lowest order i.e. dilaton, ■ • • but also by some non trivial statistical effects. 

F7I. 



We have restricted ourselves to the two extreme cases a = and (3 = 0. 
One can ask what happens when a and (3 are simultaneously non vanishing. 



An uncertainty relation similar to the one displayed in Eq. ( |3.4| ) exists but 
one does not possess a position nor a momentum representation. To deduce 
classical statistical physics from quantum statistical physics as we did , one 
has to deal with the Bargmann representation || which is more complicated. 
In this case , J will depend on both T and V. The reserves we expressed due 
to the presence of a preferred origin also apply here. 

Can a similar treatment be performed for gases in the q deformed algebra? 
The situation may prove difficult since for a free particle in this theory, the 



momentum is quantized without advocating boundary conditions [K| . 
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VI-A The Non relativistic non interacting gas in the 
a = theory 

The change of variable x = p 2 /2mkT, allows one to write J as the dimen- 
sionless integral 

/•CO 

J = 2ix~ l/2 dxe- x x 1/2 (l + 2(3mkTx)- 3 . (6.1) 

Jo 

This integral can not be performed analytically. Nevertheless a useful ap- 
proximation can be obtained which gives a good qualitative picture. Being 
the product of an increasing and a decreasing exponential, the integrand 
presents a maximum at 

-(1 + 5pmkT) + VI + UPmkT + 25/3' 2 m 2 k 2 T 2 . , 

X = 4^kf • (6 ' 2) 

Calling 4>(x) the integrand in Eq. (|3.34| ) , a saddle point approximation gives 

j = {mf 2 ^=- (6.3) 



The fact that 4>(x) is maximal at x simplifies the expression of the second 
derivative 

(j)"{x) = (1 + ApmkTx)(4:pmkTx 2 - 4x + 1). (6.4) 

In this case J is a function of the temperature but it does not depend on the 
particle number or the volume. Consequently the equation of state PV = 
NkT is still valid here but the molar specific heat at fixed volume is no 
more the constant 3R/2. Let us compute the first correction to J. A gas 
presents a non relativistic behavior at low temperature; in this limit we 
obtain, introducing the particle dependent characteristic temperature T c = 
1/ f3mk and find to the first order in T c /T 

1 T 

x = --3— (6.5) 

2 T c K J 

J = <x + ]>>(T/T C ) Z (6.6) 

with 
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a = 1.20353,ai = -10.8318, a 2 = 93.8755, 



Now 

S = Nk 



-+Wi + ln4 — 



V (2TxmkT\ 



3/2' 



fi = —kT 



In a Q + In < — 



.V \ h 2 ) 

V f2nmkT\ 3/2 



N \ h 2 J 









and the equation of state is 



'.2 ff„r c 



An adiabatic process takes the form 



V = c te NT- 3/2 e^% 



(6.7) 

(6.8) 
(6.9) 

(6.10) 
(6.11) 



while densities X t = Ni/V at equilibrium for the reaction displayed in 
Eq.(gg§ obey 



Y c Y d ( CT i t a j- b J '^ d ' 

"^C^D _ terp3/2(-a-b+c+d) I ~^ ( -^ + 73- + 7^+^p ■ 



A A A B 



(6.12) 



since each particle, having its own mass, possesses a specific critical temper- 
ature. 



VI-B The Non relativistic non interacting gas in the 
(3 — theory 

If one considers a small container, then the following approximation is accu- 
rate for the quantity J : 

J=l-a 1 aV 2/3 (6.13) 



with 



The new pressure 



4 / 3 x5/3 



-2/3 



P = NkTV' 1 (l - -a ia V 2/3 ) 



(6.14) 
(6.15) 
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is lower than the one obtained for the same gas in the unmodified algebra. 
The minus sign suggests that the pressure may become negative. However, as 
a must be tiny, this would occur for huge volumes but then the approximation 
specified in Eq. Q3.14 ) is no more valid. 



The entropy and the chemical potential become 

S = Nk 



5 T ,2/3 , I V f2TrmkT\ 3/2 ' 

- - axaV 2 ' 6 + In { — 



(jl = —kT 



N \ h 2 J 



, i V {27imkT\ ' i .,,,., 



(6.16) 



(6.17) 



while the specific heat remains the same. One sees that for small volumes 
nothing really new appears. 

VII figure captions 

• Figl The quantity 1 + g{T /T CT ) shows that the energy density due to a 
radiation of fixed pressure decreases at high temperatures in the a = 
case. 

• Fig2 The function 1 + f(x) which modifies the equation of state when 
(3 = is depicted. It tends quickly to its asymptotic value. For a system 
with fixed T, N, the pressure falls more quickly when the volume is 
increased in comparison to the undeformed case. 

• Fig3 The function fi(x) linked to the first derivative of the pressure 
does not vanish for finite volumes so that the isotherms in the (3 = 
theory do not have an extremum in the P, V diagram. 

• Fig4 The function f2(x) linked to the second derivative of the pressure 
does not vanish neither. This and the previous figure proves there is 
no phase transition. 
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